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ABSTRACT 

We present a Chandra monitoring campaign of the highly variable Seyfert galaxy 
UGC 4203 (the "Phoenix Galaxy") which revealed variations in the X-ray absorb- 
ing column density on time scales of two weeks. This is the third, clear case, after 
NGC 1365 and NGC 7582, of dramatic Nh variability on short time scales observed 
in a "changing look" source, i.e. an AGN observed in the past in both a reflection- 
dominated and a Compton-thin state. The inferred limits on the distance of the X-ray 
absorber from the center suggest that the X-ray "torus" could be one and the same 
with the broad emission line region. This scenario, first proposed for an "ad-hoc" pic- 
ture for NGC 1365, may be the common structure of the circumnuclear medium in 
AGN. 
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1 INTRODUCTION 

The structure, size and composition of the circumnuclear 
medium of Active Galactic Nuclei (AGN) is the subject of 
several studies at different wavelengths. In particular. X-ray 
absorption variability has proven to be a powerful tool in 
order to establish the physical properties of the X-ray ab- 
sorber/reflector: the almost ubiquitous absorption variabil- 
ity detected in bright absorbed AGN (Risaliti et al. 2002) 
is in itself an indication of a clumpy structure of the ab- 
sorber. The time scale of such variations provide constraints 
on the dimensions and the distance from the central black 
hole of the obscuring clouds, under the assumption of Kep- 
lerian motion. The best characterization so far of the phys- 
ical parameters of the X-ray absorber has been possible for 
the AGN in NGC 1365. This source has been monitored 
several times in the last few years with Chandra, XMM- 
Newton, Suzaku, and showed extreme variability in the X- 
rays: changes from Compton-thick to Compton-thin states 
on time scales from weeks (Risaliti et al. 2005) to days 
(Risaliti et al. 2007) to ~ 10 hours (Risaliti et al. 2009A) and 
eclipses due to Compton-thin clouds {Nh ~a few 10^^ cm~^) 
on similar time scales (Risaliti et al. 2009B). These ob- 
servations strongly support the view of an X-ray absorber 
made of clouds with Nh from 10^^ to >10'^* cm~'^, densities 
lO^'^-lO^^ cm~^ and distance from the center of hundreds of 
gravitational radii. These are all values typical of the broad 
emission line clouds (BELC). The global picture emerging 
for NGC 1365 is therefore quite different from the standard 



view of the AGN unified models (Antonucci 1993): here, the 
X-ray absorber seems not to be associated to the parsec- 
scale dusty "torus" (Krolik & Begelman 1988) obscuring the 
BELC, but to the BELC themselves. 

This scenario appears quite well justified for the single 
case of NGC 1365, whose extreme variability is so far unique 
among AGN. Two explanations of this uniqueness are pos- 
sible: (1) NGC 1365 is intrinsically dilTerent from (most of) 
the other AGN; (2) NGC 1365 is an "average" AGN with 
favourable observational conditions (e.g. a particular line of 
sight through a torus made of clouds of different Nh, as in 
Risaliti et al. 2005). 

In order to investigate this issue, the obvious strategy 
is to look for similar variations in other obscured AGN. 

The most extreme known cases of spectral variations 
in AGN are those of the "changing look" AGN (Matt 
et al. 2003), i.e. sources with past X-ray observations in 
both a Compton-thin and a reflection-dominated state. 
Members of this class, besides NGC 1365, are NGC 6300 
(Guainazzi. 2002), UGC 4203 (Guainazzi et al. 2002), 
NGC 2992 (GiUi et al. 2000), NGC 7674 (Bianchi et 
al. 2005), NGC 7582 (PiconceUi et al. 2007). In the latter 
case, a recent Suzaku campaign revealed N// variations in 
time scales of a day (Bianchi et al. 2009), making NGC 7582 
a second case of a changing look" AGN with fast Nh vari- 
ability. 

The main open question about these sources is whether 
the observed dramatic spectral changes are due to variability 
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Table 1. Spectral parameters from past X-ray observations 
Instr. Date Counts F^ V^ N 



H 



ASCA 


18/11/95 


320 


0.2"= 


2.0'' 


- 


XMM 


05/05/01 


>10* 


2.6 


2 0+0-2 


2 1+''-2 
^■1-0.2 


Suzaku 


27/05/07 


>10* 


1.4 


1.9±0.2 


q O + 0.2 
J.O_(j 2 


Swift 1 


06/10/06 


88 


1.6±0.4 


i.g'' 


^ 7+0.5 
-^■'-0.4 


Swift 2 


22/01/07 


32 


1.2±0.5 


i.g'' 


n t; + 1.5 
^■''-1.0 


Swift 3 


25/02/07 


43 


0.9±0.5 


i.g'' 


1 g+^-° 

1-^-0.6 


Swift 4 


29/04/08 


35 


2.3±0.8 


i.g'' 


04+1.2 



Results form past hard X-ray observations of UGC 4203. °: 2- 
10 keV unabsorbed flux, in units of 10~^i erg s~l cm~^; ^■. 2- 
10 keV photon index; '^■. absorbing column density, in units of 
lO'^^ cm~^; '^■. fixed. ": flux of the observed reflection-dominated 
spectrum. References: ASCA and XMM: Guainazzi et al. 2002; 
Swift: this work; Suzaku: Matt et al. 2009. 

of the intrinsic emission or of the circumnuclear absorber. 
The case of NGC 2992 is clear: multiple observations on a 
~ 20 years period strongly suggest an intrinsic fading of the 
X-ray source (Gilli et al. 2000). We note that this is the 
only case among the "changing look" sources with a column 
density in the Compton-thin state of the order of 10^^ cm~^, 
while the other sources have A'^h>10'^^ cm~^. 

Here we present the results of a Chandra monitoring 
campaign of one of the three remaining sources, UGC 4203 
(other names: MKN 1210 and "Phoenix Galaxy"). Another 
Chandra campaign on one of the last two sources, NGC 6300, 
will be discussed in a future paper. 



2 PAST OBSERVATIONS 

UGC 4203 has been observed several times in the hard X- 
rays. The two observations which led to its classification as 
a "changing look" AGN, reported in Guainazzi et al. 2002, 
were performed with ASCA in 1998 and XMM-Newton in 
2001. A further observation, has been performed simulta- 
neously with XMM-Newton by BeppoSAX. More recently, 
the source has been observed with Suzaku in 2007, and with 
Swift-XRT for four times between 2006 and 2007. In or- 
der to obtain a complete view of the variability history of 
UGC 4203, we analyzed the Swift-XRT short observations 
(about two ks each) in order to estimate the flux and, if pos- 
sible, the basic spectral parameters. We followed the stan- 
dard procedure suggested by the Swift Data Center for the 
reduction and calibration. In all cases the net 2-10 keV spec- 
tral counts are below 100, allowing only a basic spectral 
analysis. 

A summary of the spectral parameters obtained in these 
observations, assuming an absorbed power law model, is 
shown in Table 1. 



3 CHANDRA OBSERVATIONS: REDUCTION 
AND ANALYSIS 

The five 10 ks observations were performed with the ACIS- 
S instrument on-board Chandra between January 15 and 
February 06, 2008. The first three observations are almost 
consecutive, the fourth one is about two days after, and the 
fifth 17 days after the fourth. 



We reduced and analyzed the data using the CIAO 4.2 
package, and the most updated calibrations, following the 
standard procedure developed by the Chandra Data Cen- 
teijj. The background was extracted in an empty region in 
the ACIS-S field of view. The data were checked for possi- 
ble alterations due to pile-up. Despite the iiigii unabsorbed 
flux of the source (of the order of 10~^^ erg s~^ cm~^), no 
significant pile- up is expected, since the soft (< 3 keV) pho- 
tons from the AGN are absorbed by a column density of 
the order of 2-3x10^'^ cm~^. The spectral analysis has been 
performed with the Xspec 12 package (Arnaud et al. 1996). 

A visual inspection of the five spectra (Fig. lA, B) re- 
veals a soft component at E<2 keV constant within a few 
per cent, and a second, hard component at E>2 keV with 
clear, strong variations between the fifth observation and 
the others. 

We fitted the five spectra separately with a model 
consisting of the following components: 

1. a soft spectrum produced by coUisionally-ionized diffuse 
gas, with two different teniperatures and densities (modeled 
through the APEC modeO in XSPEC), and a power law, 
representing a possible scattering component from the 
central AGN. Our aim is to obtain a statistically acceptable 
representation of the soft emission, in order to fix it in 
the analysis of the variability of the hard component. 
An alternative, simpler approach would be an analysis 
of the spectrum at energies E>3 keV, neglecting the soft 
part of the spectrum. This would however prevent us 
from rejecting models reproducing the hard spectrum, 
but whose extrapolation at low energies over-predicts 
the observed counts. The fitting analysis shows that 
all the three components are needed to obtain a good 
representation of the observed data. The thermal compo- 
nent has a temperature kT=0.7±0.2 keV and the power 
law slope and normalizations are F^ = 2.05 ± 0.15 and 
Ns = 2.0^^-3 xlO"'^ keV s"^ cm'^ keV"^ (corresponding 
to a flux of a few per cent of that of the intrinsic AGN 
continuum, see below). These are all typical parameters 
for the soft emission of local AGN. The metallicities of the 
photoionized component span a large range (from 0.5 to 
20 solar) and are poorly constrained. This is analogous to 
what is obtained in global fits of the soft diffuse emission 
of other AGN. A spatially resolved analysis of these cases, 
when possible (i.e. in lower-redshift sources) show that 
the total soft emission is the result of the contributions of 
different zones in the galactic bulge with a large range of 
metallicities and densities (e.g. Baldi et al. 2006, Wang et 
al. 2009). Since we were interested in the hard component, 
we froze the soft components to their best fit parameters in 
the subsequent analysis. 

2. An intrinsic AGN emission consisting of an absorbed 
power law. The absorption is modeled with a single pho- 
toelectric component, and then with a double component, 
one of which partially covering the source (see below for 
details). 

3. A cold reflection component, consisting of a continuum 
estimated with the XSPEC PEXRAV model (Magdziarz & 



1 http://cxc.harvard.edu/ciao/index.html 
^ http://cxc.harvard.edu/atomdb/ 
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Table 2. Results from the spectra analysis of Chandra data 
Group r N"^ R'' F" xVd-o.f. 







Independent 


fits 






OBS 1-3 


1.83«-i? 


32+^ 


,7-1-1.7 
^•'-0.8 


2.2 


176/215 


DBS 4 


0.6412:^3 


i4.7+_l:l 


<7 


1.1 


62/64 


OBS 5 


77-1-0.66 

^■"-0.68 


211- 


<8 


1.1 


51/49 




Constant reflection, 


1.7<r<2.2 




OBS 1-3 


1.7-2.2 


29. 61^-* 


1 c-l-i.o 
^■0-1.0 


1.8 


306/326 


OBS 4 


1.7-2.2 


25.5+^-^ 


2.0t!-^ 


1.4 


306/326 


OBS 5 


1.7-2.2 


37.6ltj 


i-9t^:^ 


1.4 


306/326 



■^ Column density in units of lO'^^ cm~^. '': Ratio between the 
intrinsic continuum assumed in the reflection component and the 
observed continuum. '^: unabsorbed 2-10 keV flux, in units of 
10~^^ erg cm~^ s~^. 

Zdziarski 1995) and an iron emission line. 

Since the spectral fits confirmed that no significant vari- 
ation is present among the first three almost consecutive ob- 
servations, we repeated the analysis linking the parameters 
in these intervals. Hereinafter, we use the names OBS 1-3, 
OBS 4, and OBS 5 to refer to the three groups considered 
in the analysis. 

The results from the spectral fitting are reported in Ta- 
ble 2 (top). Even if the models are acceptable from a statis- 
tical point of view {xt ~ 1 for all the spectra), the present 
analysis is not satisfactory for two reasons: 

1) the best fit continuum in the fourth and fifth interval 
is extremely flat (r~0.7). Such a flat photon index has not 
been observed thus far in any non-reflection dominated AGN 
(see, e.g. the distribution of F in nearby AGN in Bianchi et 
al. 2009A, or Dadina et al. 2008), and is therefore unlikely 
to be real. Moreover, all the previous spectra of UGC 4203 
with high statistics (Table 1), and the combined Chandra 
spectra of the first three interval, also having a much higher 
signal-to-noise, are well fitted with a fairly "normal" photon 
index, in the range 1.7-2.2. 

2) The spectral parameters in OBS 4 and OBS 5 are poorly 
constrained, due to the high degeneracy between photon in- 
dex, column density and reflection component. This is evi- 
dent from the large errors reported in Table 2. However, the 
F-N// contour plots (Fig. 2), show that a constant solution 
for all the intervals is ruled out at high statistical signifl- 
cance. 

In order to overcome these limitations and better con- 
strain the spectral parameters, we repeated the spectral 
analysis fltting all the interval simultaneously, and requir- 
ing a) a constant reflection component, and b) a continuum 
photon index in the range 1.7-2.2. Not surprisingly (given 
the large contours in Fig. 2) the overall flt is completely 
satisfactory (x^=306/324 d.o.f). The results are listed in 
Table 2 (bottom) and shown in Fig. 3. Both variations in 
intrinsic flux and in absorbing column density are necessary 
to explain the observed spectral variations. With these pre- 
scriptions, treating the three column densities as indepen- 
dent, the probability of having no column density variation is 
P(ANif=0)<10~^). Speciflcally, the measured column den- 
sities in OBS 1-3 and OBS 4 are consistent, while the column 
density from the flfth observation, N_ff (OBS 5), is different 
from N_ff(OBS 1-3) at a 99.5% confidence level, and from 




Energy (keV) 

Figure 1. Chandra 0.5-10 keV spectra of UGC 4203. Note that 
the soft (E<2 keV) and hard (E>7keV) emission is constant 
within a few percent in all the flve observations, while signifl- 
cant variations occurred in the 2-6 keV range between the flrst 
four observations and the flfth one (in black). 
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Figure 2. Photon indcx-N^f 1- and 2-cr contours for OBS 1-3 
(red), OBS 4 (green) and OBS 5 (blue). 



Nj:/(OBS 4) at a 99.98% confidence level. The latter is the 
statistically most significant difference. The contour plot for 
Nj:/(OBS 4) and Nh(OBS 5) are shown in Fig. 4. 

Finally, repeating the fit assuming constant col- 
umn density between OBS 1-3 and OBS 4, we 
obtain Nh(OBS 1-A)^29.2t{ixl0'^^ cm'^ and 
N//(0BS 5)=37.7+4gXlO^^ cm~^ The probability 
of no column density variation is obtained imposing 
Nj:/(OBS 5)=Nff(OBS 1-4). P(ANh=0)=3x10-^ 

We conclude that the column density variability be- 
tween OBS 5 and the previous intervals is highly significant. 
Changes in intrinsic flux are also present, and the contin- 
uum spectral slope may be also variable (within a physi- 
cally acceptable range), but these components alone cannot 
reproduce the observations without a variable absorber. 
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Figure 3. 3-10 keV spectrum, best fit model and x^ residuals for 
the five Chandra observations. 
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Figure 4. Column density contours (1, 2 and 3 a) from tlie best 
fit models to OBS 4 and OBS 5. The equal column density line 
is also plotted, for ease of comparison. 



4 THE STRUCTURE OF THE X-RAY 
ABSORBER 

The results of our analysis show a column density varia- 
tion in a time scale AT between 3 days (the time between 
the first and the fourth observation, during which no vari- 
ation has been found) and 17 days. This result presents a 
strong resemblance to the column density variations found 
in the analogous Chandra campaign on NGC 1365 (Risaliti 
et al. 2007). Therefore, in the following we apply the same 
analysis as for NGC 1365. 

The variability time scale can be easily translated 
into an estimate of the distance of the absorbing cloud, 
assuming a rotation around the central black hole (BH) 
with Keplerian velocity Vk- This estimate is subject to 
large uncertainties, connected to the values of the BH mass. 



and of the cloud velocity, which are discussed below. 
Black hole mass. An estimate of the BH mass Ms// 
can be obtained from the known relations between M-bh 
and spectral lines, bulge luminosity, and bulge velocity 
dispersion. In the case of UGC 4203, discrepant estimates 
are found in the literature: 

1) MflH-o" relation. The relatively low bulge velocity disper- 
sion (ct=77 km s~^, Garcia- Rissmann et al. 2005) suggests 
log(MsH/M0)~6.5 (Bian & Gu 2007). The value obtained 
from the velocity dispersion imply a super-Eddington 
luminosity, if a standard X-ray to bolometric correction 
is adopted, based on known Quasar Spectral Energy 
Distributions (e.g. Risaliti & Elvis 2004), and the relation 
between optical to X-ray ratio and UV luminosity (Young 
et al. 2009). We consider this scenario highly unlikely, 
mainly because of the rather constant intrinsic X-ray flux 
(within a factor ~2 of the one measured with Chandra, 
Table 1 and 2) over several years, which would imply a 
persistent super- or near-Eddington luminosity for this long 
period. On the other hand, the estimate of Ms/r based on 
velocity dispersions may be strongly biased in this case by 
the presence of an unusually strong starburst in the bulge 
of UGC 4203 (Mazzalay & Rodriguez- Ardila 2007). 

2) Ms/r-bulge luminosity relation. We estimated M-bh 
using the Marconi & Hunt (2003) relation between H-band 
bulge magnitudes and BH masses, using the 2MASS H 
magnitude in the inner 4 kpc of the galaxy, and the 
ratio between AGN and stellar contribution in H-band 
estimated by Mazzalay & Rodriguez- Ardila (2007). We 
obtain log(MflH/M0)~7.7. 

3) Virial mass estimates. The reverberation mapping- 
calibrated relation of Vesteergard et al. (2006) suggests 
log(MB/i/Mo)~7.9 (Zhang et al. 2008). In this estimate 
the width of the broad H/3 line is measured in the polarized 
spectrum, while the intrinsic luminosity at 5100 A is 
estimated from the [O HI] A5007 A hne. - The [O III]-based 
luminosity is confirmed by an analogous estimate based 
on the intrinsic X-ray luminosity, following the same 
prescriptions as mentioned above. The measurement of 
FWHM(H/3) is also quite solid, being based on a rather 
high quality spectrum (Tran 1995). 

Concluding, despite the large uncertainties in each 
of the methods summarized above, we assume MB/r~5- 
7x10^ niQ as the most likely mass estimate. 

Cloud velocity, the linear dimension of the X-ray 
source cannot be lower than a few gravitational radii, 
_RG=GM_B/r/c^. Adopting Ds'^5 Kg as a lower limit, the 
velocity of the eclipsing cloud must therefore satisfy the 
condition: V/s-xAT>5 Kg- Adopting the mass estimate 
discussed above, this translates into Vk>10^ km s~^ 
ATio, where ATio is the time interval in units of 10 days. 
This value is clearly in the range of the Broad Emission 
Line Clouds (BELCs) and well above that typical of the 
standard, parsec-scale obscuring torus. 



4.1 Physical properties 

Based on the above estimates, and assuming Keplerian mo- 
tion of the obscuring cloud, we obtain a distance from the 
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central black hole R~l-2x 10^^ cm. The cloud linear size Dc 
must be of the order of, or a few times larger, that that of 
the X-ray source, i.e. Dc~l-5xlO^* cm. The measured col- 
umn density Nh^W^^ cm~^ imply a density of the order of 
lO'* cm-^ 

We conclude that the X-ray absorber has all the typical 
physical parameters of a Broad Line Region cloud. There- 
fore, absorption in X-rays and the broad optical/UV emis- 
sion lines are most likely due to the same physical compo- 
nent. We also note that UGC 4203 is optically classified as 
a type 2 source, implying the presence of dust along the line 
of sight. This dust cannot be associated to the Broad Line 
Region, which is well inside the dust sublimation radius, 
and must therefore be located in a farther absorber. The 
nature of this second absorber may be a galactic dust lane, 
and/or a "classical" parsec-scale torus. In either case, a con- 
tribution to the observed X-ray absorption by this second 
component is expected. Within the limited statistics of our 
Chandra observations it is not possible to disentangle these 
two components in the spectral fits. The contribution to the 
observed column density may be as low as a few 10^^ cm~^, 
or as high as the lowest Nh measured in past observations, 
i.e. ~2xl02^ cm-2 (Table 1). 

UGC 4203 is one more "Changing look" AGN where 
the extreme X-ray variability is proven to be due to absorp- 
tion by BLR clouds. The same conclusion can be reached 
for a few more sources where a similar variability has been 
found (NGC 4151, Puccetti et al. 2007, NGC 4388, Elvis et 
al. 2004, NGC 7582, Bianchi et al. 2009B, Mrk 766, Risaliti 
et al. 2009). On the other hand, it is likely that other bright 
sources with multiple X-ray observations do not show short- 
term absorption variability (even if no quantitative analysis 
has been done so far), and that other sources of X-ray ab- 
sorption, quite far from the central source, may be present 
(e.g., as discussed above, galactic dust lanes and/or parsec- 
scale tori). Therefore, since the physical link between BLR 
clouds and X-ray absorber has been now established for a 
good number of sources, our next step will be the complete 
analysis of a representative sample, in order to quantita- 
tively estimate the relevance of short term variable X-ray 
absorption in local AGN. 
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